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ABSTRACT 


Evolution of thin viscoelastic liquid films has been studied using Oldroyd-B model for stress- 
strain rate constitutive relationship and the simplifying the Goverening differential equation 
under lubrication approximation. 2-D Numerical simulations have been performed of the 
nondimensionalized equations. Presence of a polymer in a Newtonian solvent makes the solution 
viscoelastic. Instead of being linear, long polymer chains try to bring each other to randomize the 
system and hence increase the entropy, which shows the elastic behavior of viscoelastic fluids. It 
was found that the length scale of instability is not affected by the inclusion of elasticity in the 
liquid, which is due to the presence of polymer in the solution. But time scale of instability is 
much affected and is very sensitive to the elasticity of fluid especially in Maxwell model when 
only polymer is present in the solution. As the elasticity of solution is increased there is a 
decrease in the time scale, which shows the destabilizing effect of elasticity. As the film surface 
touches the substrate and starts forming a hole, a rim around the hole fonns and rim height 
increases linearly with the hole diameter. Effect of bending moment in the Newtonian fluid is 
also considered. Bending moment, like the surface tension, has a stabilizing effect. 3D 
simulations are also carried to show how the pattern changes with the height. 
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CHAPTER 1 


INTRODUCTION 


I. 1 Thin Films 

In thin films the thermodynamic and transport properties depend on the film thickness, 

which is less than the effective range of intermolecular interactions (usually < 100 nm). 

Thin films have much importance not only in scientific interests but also in the industries. 

Stability, morphology, adhesion failure, and defect formation in thin films are of main 

scientific interest in the areas like coating, painting, adhesives, floatation, multilayer 

adsorption etc. Thin films are widely used, for example, to provide insulating layers 

between conductors, diffusion barriers, and hardness coatings for scratch and wear 

) 

resistance. Excess intermolecular interaction, which is the characteristic of thin films, is 
the main cause of the thin film instability. Excess molecular interaction is the sum of 
many possible components, the most important of which are Lifshitz-van dcr Waals 
forces (LW), electrostatic (EL) interactions, polar “acid-base” (AB) interactions and the 
extremely short range Bom Repulsion (BR). Depending on the combination of these 
force fields, the thin films are characterized to four types of films namely. Type I, Type 

II, Type III, Type IV. The force fields in Type I, Type II, Type III and Type IV films are 
of purely attractive, long range attraction-short range repulsion, purely repulsive and long 
range repulsion-short range attraction respectively. 



1.2 Literature Review 


Instability of liquid thin films has been studied extensively and systematically for 
more than 40 years. Many possible pathways have been proposed and demonstrated for 
the transition from the initially flat and smooth film to the final state consisting of 
droplets. The free surface of an initially uniform thin film becomes unstable and deforms 
spontaneously to engender a microstructure when the second derivative of the excess 
intermolecular free energy (per unit area) with respect to the local film thickness is 
negative, viz., <0[l-5]. Investigation of opening of holes in thin films, 

deposited onto silicon wafer coated with a polydimethylsiloxane monolayer, shows 
strong deviation from the behavior of pure Newtonian liquid. The highly asymmetric rim 
grows which varies with the diameter, linearly [11]. Polymer melts, which are Non- 
Newtonian in nature, do not obey linear relationship between shear stress and strain rate. 
Many models have been proposed to model the constitutive relationship between stress 
and strain of such fluids, namely. Maxwell fluid model, Voigt solid model and Oldroyd-B 
model. Oldroyd-B model is most general and linear viscoelasticity is considered here. 
The length scale of instability, in Non-Newtonian fluids, is unaffected of viscoelastic 
behavior of the fluid and it remains same as in the Newtonian fluids. Viscoelasticity does 
not determine the competition between the capillarity and intermolecular potential that 
set the scale for a viscous fluid while the time scale is decreased. 

S. A. Safi-an and J. Klein [12] state that the surface of liquid thin films can be unstable 
due to thinning van der Waals interactions, leading to the formation of holes in the 
initially uniform film. These instabilities can be greatly retarded in viscoelastic materials 
and completely inhibited in the elastic materials. 
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A. Sharma et al [6] have shown the formation and growth of holes in apolar non-slipping 
Newtonian thin films subject to long-range Lifshitz-van der Waals forces are investigated 
based directly on numerical solutions of the thin-film. The present study is an extension 
of this work, where the fluid is considered to be Non-Newtonian in nature. 

1.3 Outline of the present study ^ 

The chapter 2 deals with the instability in viscoelastic thin films on a 
homogeneous substrate. Figure 2.1 shows the viscoelastic fluid on a substrate, bounded 
by a bounding media like air. Some examples of viscoelastic fluids include toothpaste, 
soft serve ice cream, ketchup, melted chocolate and tar. We start with the theory of 
viscoelastic fluids, Oldroyd B model, and then constitutive equation is derived. Then for 
this 2-D system, Navier Stokes equation is applied, with the appropriate boundary 
conditions, which is scaled using the long wave approximation. The instability of the 
liquid film is studied, resulting, evolution of the free surface as a function of time and is 
compared with the linear theory. The length scale and the time scale of instability are, also 
compared with the Newtonian case. 

Chapter 3 deals with the Newtonian fluid on a solid substrate bounded by the cell 
membrane. The effect of bending moment on the cell membrane is studied and its effect 
on the length-scale and time-scale with the change in intensity of bending moment has 
been shown. 3-D simulation performed also show the variation in patterns with variation 
of the mean height of the film. 
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CHAPTER 2 


INSTABILITY IN VISCOELASTIC THIN FILMS 

2.1 Theory 

Oldroyd B model is one of the many models, to understand the behavior of viscoelastic 
fluids, which has both the viscous property of liquid and elastic property of solid. In the 
Newtonian solvent, which has low molecular weight, the polymer is present which makes 
the solution to be viscoelastic, polymers are long, randomly arranged chains so that it can 
be considered as two balls connected by an elastic spring. To understand the elastic 
property of a polymer in the solution, we can think of long polymer chains in the solution 
to be randomly arranged, due to entropy constraints and have a tendency to pull back to 
make the process more random, and hence, increase the entropy. While doing so, the 
viscous behavior of the solution and polymer makes the process slow. 
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^G, 

Figure 2.1; shows the mechanistic equivalent of Oldroyd B model, Spring shows the 
elastic property of polymer with spring constant Gi. 


2.2 Constitutive equation 

In a Newtonian liquid, stress t is proportional to the strain rate ^ , i.e. to the gradient of 
the velocity filed^^ > but in a linear elastic medium the stress is proportional 

to the strain ^ , but not the strain rate. The total displacement ^ in the Oldroyd B model is 
the sum of displacement due to polymer and Newtonian solvent and respectively. 

Stress in spring is proportional to strain while in dashpot, it is proportional to strain rate. 
Hence total stress can be written as the summation of these two forces. 




\dt j 


,(2.1) 


Similarly, stress in the polymeric dashpot can be written as. 




p_ 

V J 


( 2 . 2 ) 


The total displacement is the sum of these two displacements. 


( 2 . 3 ) 
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the constitutive equation for the Oldroyd B model is given by 


dt 




A 7 , 




dt 


V ^ J 


dy 


(2.4) 


Where, x is the direction parallel to the film surface and y is the normal direction to the 
film surface, X. is the relaxation time (the time at which forces relax to its 37 percent)[9]. 


5F , 


dy 


is strain rate (i^) and r] is the sum of viscosities of Newtonian solution and 


polymeric melts (7 = 7 ^. + 7 ) and 7 ^ is the ratio of solvent viscosity to the total viscosity. 


2.3 System of Study 

Here we consider the thin viscoelastic liquid film (<100 nm.) --polymer mixed with 
Newtonian solution, - which overall makes the solution viscoelastic, bounded by a solid 
support (perfectly homogeneous) from below and a bounding media (air for example) 
from above, is the subject of the current study. The length scale parallel to the film 
surface is much larger than the length scale perpendicular to the film surface and hence 
the lubrication approximation is valid. Such a system can be better visualized with the 
help of Figure 2.1 which is a schematic diagram of a supported thin viscoelastic liquid 
film. Here we are considering only two dimensions one parallel to substrate and one 
normal to the substrate. 
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Bounding medium 



Substrate 

Figure 2.2 Schematic diagram of a 2-D thin film on substrate. 


2.4 Free Energy of Thin Films 

The equilibrium and dynamics of thin fluid films, like that of any other physical system, 
are related to the changes in the free energy rather than the free energy itself. Thus to 
understand the' changes in free energy is of paramount importance to describe the 
mechanics of thin films. The contribution to the free energy in a system is because of 
interfacial tension, but if the film is less than 100 nanometer, then interrholecular forces 
also play an important role. Obviously the free energy change also depends on the film 
thickness. The disjoining pressure and excess free energy are obtained by pair wise 
summation of interactions among molecules of the thin film and among molecules of the 
film and the solid. The use of Lennard-Jones potential for intermolecular interactions 
gives the free energy per unit area AG and the conjoining pressure as the sum of their 
respective components due to the long-range Lifshitz-vander Waals attraction and the 
extremely short-range Bom repulsion (5, 6, 8). Under conditions of small slope, these are 




(2.5) 


And 



85 


( 2 . 6 ) 
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is the cut-off distance i.e. at which (j) vanishes. And we know \h 2 Li,A = -\2KdlS ‘''^ , 
replacing these values in the above equation and using the non-dimensional parameter in 
equation (2.1 1) one can get the non dimensional equation for disjoining pressure as, 


n° = 


d-H 1 


dX- 3//' 


1 - 


r I \ 






(2.7) 


This fomi of potential belong to type I films where free energy change AG and its second 


d-G 

dH^ 


derivative, is negative for all film thickness. 



Figure 2.3. Qualitative variation of change in free energy AG and second derivative of 


3^(j 

change in free energy , with the film thickness for a type I thin film. 

dH 


2.5 Conservation laws 




Here we will state the hydrodynamic equations of viscoelastic fluid, and will start with 
the conservation laws. 
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2.5.1 Mass conservation 


Assuming the density of fluid to be constant the equation of mass conservation reduces to 

V.u = 0 (2.8) 

Where V is the velocity vector (Vx, Vy). 


2.5.2 Momentum conservation 

p — = -Vn + V.T (2.9) 

dt 

Here, — = 9, + u.V is the total (or material) derivative. And fl is the total pressure which 
dt 

comprises the pressure inside the film due to interfacial tension and the excess pressure, 
(f ) , which is negative of the disjoining pressure, and is related to the excess free energy 
per unit area of film as 




aAG 

dh 


( 2 . 10 ) 


The body force due to gravity is neglected compared to the body force due to excess 
intermolecular interactions. 


2.5.3 Boundary Conditions 

At the free film surface tangential component of stress is assumed to be zero and 
at the film-substrate interface, no slip boundary condition is assumed. Also the wall is 
assumed to be impermeable to the fluid, hence normal component of velocity at the 
substrate surface is assumed to be zero. 
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2.6 Scaling 

The following parameters have been used for non-dimensionalizing the governing 
equations. 

f nh \ 

H = h/h, T = !/ ( 2 . 11 ) 

I r J 
(y^ 

V = u/ — T = T I 

\n) 

• Hq is the mean film thickness 

• 77 is the sum of both the viscosities. 

, • 7 is the interfacial tension of liquid 

Substituting these dimensional variables in terms of non-dimensional parameter in the 
above constitutive equation 2.4, we get the constitutive relation as 



■ + 'r.yy =w,Vr- 


d fdK) dV. 


dT [dr J dr 


( 2 . 12 ) 


Tyy is same as r in equation ( 2 . 1 1 ) and w. is the weissenberg constant, defined as 



(2.13) 


For very thin film the length scale of film thickness \ is very small than the length scale 

L parallel to the substrate surface. Hence, under the lubrication approximation the inertial 
term is not important and only the viscous term is comparable to the pressure term. 
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Navier Stokes momentum equation, in the x direction, under the lubrication 
approximation reduces to, 


an _ dT^. 

dX dY 


(2.14) 


and. 


m 

dY 


= 0 


(2.15) 


From the equation (2.15) we can say that, n is function of X only. And hence equation 
(2.14) can be easily integrated. Integrating the equation (2.14) with respect to Y we get 


stress as 




(2.16) 


Substituting the value of stress, obtained from equation (2.16), in the equation (2.12) we 
get 


w, 


' dT 


— (X))l+— (JF)) = — 

^ dx^ ^ dY 


dVy 


(2.17) 


Further integrating the above equation with respect to Y we get 


w. 


' dT 
(2.18) 


an 


/ v2 


dX 


■-Y-H{X,T) 


dnf Y- 


dX 




-Y-H {X,T) 


J 


dVy 

dT 


(0) 


Renormalizing the n we get 

n=e^n 

Here s is the ratio between mean height of film and length scale. 


(2.19) 


Re-ordering the equation (2.18) we get. 
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L—y-h(x,t) + ^ 1—y-h(xj) 


‘ dT I dX {2 

( 0 ) 

oT 


dX I 2 


( 2 . 20 ) 


Further redefining the variable 


^ ^ ^ 


( 2 . 21 ) 


Putting these redefined values in equation (2.20) 




d'f L aA' 2 

- , 0 ) 

JJX’, -^ + v'°' 

oT 


dX 2 


( 2 . 22 ) 


For the 2-D flow, the kinematics equation can be written as 

— + — I (f) = 0 

dT dX ^ ■' 


(2.23) 


dH a _ y. , . ^ 

— + — /( 2 r,r) = o 
ar dx ^ ’ 


(2.24) 


Mere / is the total flux passing per unit width of the film. Putting the value of Y = H in 
the equation (2.22) we get the equation in interfacial velocity 


- ^Yv./ , ,0, a fan® i an° 

11 X’, —z:;r^^x i = — H' H' 

dT 2 dT dX 2 dX 


(2.25) 


Here is the interfacial velocity in the x direction. 
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Multiplying the equation (2.22) by the dY and integrating from 0 to H(X, T) and also 
using the Leibnitz theorem, the equation (2.22) simplifies to 


— r:= + 

ol dT 


w, d 
'T~dT 


an° 

dX 






2 dT dX 


3 dX 


(2.26) 


Let 


— = w • and f = — and hence — 
3 3 dT 


using these renormalized parameter in 


equation (2. 24), (2. 25), (2. 26) we get 3 sets of coupled partial differential equations, which 
governs the instability and spatio-temporal evolution of a non-slipping viscoelastic thin 
liquid film subject to excess intermolecular interaction, as 


||+3A/(x,r) = 0 (2.27) 


, di{x,T) , a 

^ j ^ + 3 / (X,T) = -w-^ 

dT dT oT 


an° 

dX 




3wa//an°„, an° j 
2 dT dX dX 


(2.28) 


3v,. , w a 

TJW , = 

dT • 2 dT 


an” 

dX 


H- 


1 an” 

— n 

2 dX 


Where, 


an” 

dX 


_a_ ^ 

dX dX^ 


1 

3H^ 


f 

V 




U 


n6 > 

^ J 


I 


(2.29) 


(2.30) 


The parameter used to non-dimensionalised the equation used are as follows. 
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// = A, A' = 

K 


X 


and 



r 


1/2“ 

/7,(/7„/t/o) 

_ z 


- 







<6 = ,^ V 

[6(-5“)<i;/*.'] 

the liquations (2.27), (2.28) and (2.29) arc the three coupled partial differential in H , 
I and V ^ / . Here / -1^/ where is the separation distance, where excess free energy 

an" 

is minimum. is the derivative of disjoining pressure, the first part, second 

dX 

derivative in height, shows the interfacial pressure term while the other part, which is 
intennolecular forces are made up of two temis, first is LW component of Vander wall 
forces and second is due to Bom Repulsion. The viscous effect in thin films merely 
retards the growth of instability. The surface tension effect has stabilizing effect because 
of the in-plane curvature at the film surface. 


• w is measure of elasticity, higher the value of w , higher the elasticity. 

• T]^ is the ratio of solvent viscosity to the total viscosity, if this tend to zero means 

we get simply the polymer solution and if it tends to 1 means we get the 
Newtonian solution. 


2.7 Linear Stability Analysis 

The initial growth of the instability in thin films can be described by linearizing the above 

governing equation, wherever the amplitude of the initial disturbance is much smaller 
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tnan tne mean tilm thickness. Height is being linearized around H=1 and other parameter 
like / and v are linearized around zero (initially there is no interfacial velocity and mass 
flux) as: 


H = 1 + 1 ^ exp{iKX + sT) 

(2.31 a) 

I = JexpiiKX+sT) 

(2.31 b) 

V = V exp{iKX + sT) 

(2.31 c) 


In the above set of equations ^,K,T are the non-dimensional amplitude (i.e. ratio of 
initial disturbance amplitude and mean film thickness), wave number and growth 
coefficient of the initial disturbance, respectively. Linearizing the governing partial 
differential equation using these parameters we get a relation between non-dimensional 
growth coefficient and non-dimensional wave number. 

7,W+{l-w(/f^-A:')}5-(X--if") = 0 (2.32) 

Also, the disturbance in velocity and flux can be obtained in terms of disturbance in 
height and hence fiux and velocity can be written as: 

H = \ + <^expiiKX + sT) (2.33 a) 

1 = ^exp{i{KX + 7r/2) + sT} (2.33 b) 

3K 


eKp\i(KX + 7r/2)-¥sT] (2.33 c) 

2 (1 -I- ;7^w5') 

The above equation which is a quadratic equation in‘s’, can be solved to get the linear 
dispersion relation, as a function of non-dimensional wave number. The critical wave 
number, at which growth coefficient is zero, comes out to be = 1 i.e. only the modes 
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with non-dimensional wave number less than unity can grow, and the non-dimensional 
length scale of the fastest growing linear mode, where growth coefficient is maximum, 
comes out to be I'Jln . , 

The maximum value of growth coefficient is: 


.5 


tnax 


-4 + w + -sj\6 + 8(2/;,. - l)vi' + vv' 


(2.34) 


Further more, the time of mpture (i.e. the first time at which film surface touches the 
solid substrate) from the linear stability analysis comes out to be: 


T =■ 


■log. 


^1^ 




(2.35) 


2.8 Numerical Method 

The following (2-D) set of Non-dimensional equation, governs the stability and spatio- 
temporal evolution of the thin film subject to the excess intermolecular interactions. 


f-3A;(^,r) = 0 (2.36) 


, di(x,T) , a 

^ 2r}^wv^,^ + 3I(X,T) = -w^ 


dX 




3wdHdn\^. 5nV,3 

-I- =r H~ H 

2 dT dX dX 


(2.37) 


dv 


7rW- 


X,1 


dT 


+ ^A'./ 


w d 


an" 

a^r 




2 dX 


(2.38) 
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3 

These are coupled partial differential equations with — appearing in the right as well as 

dT 

in the left side of the equation; i.e. the above equation is not explicit in time. Because 

these are not explicit differential equations, hence all — ^ terms are collected in one side, 

dT 


and is replaced by ^he last two equations, to make it explicit in space. 

This is fourth order in space and first order in time. We convert the above Partial 
differential equation into a set of ‘ W ’ Ordinary differential equation, using the finite 
difference technique. These equations are discretized in space using central differential 
method with half node interpolation. Initial condition on height can be taken to be either 
sinusoidal or random perturbation; while for interfacial velocity and flux are assumed to 
be zero in the initial stages. Collecting the time derivative terms in one side, the equation 
converts to: 


S + .l — /(AC7’)={) 

dT dX ^ ’ 


di 9w di an“ an° „3 


~^(a) = -'il , 

dT^ '^ ^ dX 2 dX dX 

d . \ 1 5n° 2 

ar^ 2 dx 




Where, 


1 r „3 

p, = 3rj^wI + w—H 


P2=t?,m>v^ + 


w dU° 
2 dX 




dX 


■W 


( 2 .. 19 ) 

( 2 . 40 ) 

( 2 . 41 ) 


( 2 . 42 ) 

( 2 . 43 ) 
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The above set of partial differential equations was solved with periodic boundary 
condition over an interval sufficiently large compare to 2^/2;r . The initial condition was 
small amplitude = 0.01) sinusoidal or random perturbation. NAG library subroutine 

D02EJF, which employs the Gear’s algorithm and is most suited for solving coupled, stiff 
ODES, was used to solve these set of ODES. 700 grid points in 2A scale is found to be 
sufficient for grid conversion test. 
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2.9 Results and Discussion 


Like the Newtonian fluid in viscoelastic fluids also the three distinct regimes of the hole 
growth are identified. 

1 . A short, unsteady phase in which there is faster change in velocity as well as in 

contact angle. j 

2. Quasi-steady phase in which changes are very slow. 

3. Finally the rims of a hole starts impeding the neighboring hole rim and form a 
bigger equilibrium drop. 

2.9.1 Random perturbation 

Figure 2.4. Shows the evolution of instability starting from nearly uniform film 
thickness (with a random perturbation of(£ = 0.0 1) ) until the formation of a dry drop, 

i.e. the first time when film surface touches the substrate surface. A length scale for 
the simulation was chosen of 3A , with the initial film randomly perturbed , and in the 
lust phase of evolution 3 waves started forming. Hence the length scale predicted by 
the simulation is in compliance with the length scale predicted by the linear theory. 
Length scale predicted by linear theory is the same as we obtained for the Newtonian 
case. But the time scale of instability (will be shown later) is lesser than the 
Newtonian case. The instability occurs because derivative of conjoining pressure 
with respect to height is negative, which shows higher pressure at lower thickness and 
vice-versa. In the beginning, if a random perturbation is given, then perturbation 
rearranges itself and in rearranging its amplitude decreases initially and after some 
time it starts growing and then ruptures, leading to rim interaction and finally droplet 
formation takes place. 
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Ingure 2.4. Various stages of evolution of instability in an initially uniform film thickness 
of (/?o = 6.85) nm. And weissenberg number(w = 4.0) , on a length scale of3A . Curves 1- 
5 correspond to Tr = 0, 12.282, 15.906, 17.307, 18.310. 

Small slope approximation is valid as long as the contact angle between liquid and solid 
at all stages of the hole growth is small. Visually it gives the impression that at the each 
stage of the hole growth, the contact angle is high but this is because we are taking larger 
value of X axis domain on a smaller scale in plotting, compare to the thickness scale. 
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2.9.2 Sinusoidal perturbation 

Figure 2.5 shows the full morphology of instability in thin viscoelastic liquid films under 
the sinusoidal perturbation, for 7 , =0.5, w = 4.0. The wavelength of the sinusoidal 

perturbation is taken as the dominant wavelength and simulation is carried over the 
length scale of two-wave length. The evolution comes out to be symmetric i.e. both the 
waves start growing with the same rate, touches the surface, hole growth and its evolution 
all are occurring at the same time. If the simulation is carried out for the wavelength less 
than the critical wave length (predicted by linear theory) i.e. at which growth coefficient 
is negative, it was found that perturbations does not evaluate and after some time all the 
perturbations die out, which matches well with the linear theory. A rim forms rapidly 
around the expanding hole (curve 4) and grows with the expansion of the hole (curve 5 
and 6 in Figure 2.5). As the rim starts to grow, it feels the effect of neighboring rims also 
(curve 7) which causes hindered phase of expansion. And finally merger of rims takes 
place and they form a bigger drop. Film evolution until rupture is slow, after rupture, the 
time taken to rim formation, coalition and finally forming a bigger drop is very less. 
Faster evolution after rupture is due to the surface tension which although initially was 
opposing the instability but later supports the formation of the drop, which corresponds to 
another more stable (meta-stable) film configuration. 
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Figure 2.5. Full morphology starting with a sinusoidal perturbation in a viscoelastic film 
of initial thickness = 6.85 and [ t ]^ = 0.5,w = 4.0) over a length scale of 2A , curve 1-8 
corresponds to ta = 0, 5.953, 8.666, 9.019, 9.021, 9.026, 9.037, 9.0482. 
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Figure 2.6 shows the variation of maximum growth coefficient as a function of 
weissenberg number with 77 ^ as the parameter. When 77^=0 simply we get Maxwell fluid. 

The variation of maximum growth coefficient as a function of weissenberg number for 
Maxwell case comes out to be, 

Ifw<4 (2.51) 

4— w 

And for Newtonian case ( 77 ^ = 1 .0) it is independent of weissenberg number. 

For Newtonian fluid (2.52) 



Figure 2.6.Variation of maximum growth coefficient as a function of weissenberg 
number, 77 ^ as the parameter. 
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Some important results as predicted by linear theory are given in the following paragraph. 

The dominant wavelength is proportional to (-5 Vhich means that for thicker 

films and less non-wettable substrate would respond with longer dominant waves during 
the evolution of instability. The number density of holes is proportional to 

/zq"’ ) which means that thinner films and more non-wettable films have higher hole 


density. The minimum time of mpture is proportional to 



which means that 


thinner films and higher non-wettable films are prone to mpture faster. The cut-off 
distance and dominant wavelength predicted by linear theory are also in compliance with 
the nonlinear numerical solutions. 

Figure 2.7 shows the amplitude as a function of time from linear theory and non-linear 


simulation. 




Figure 2.7. The deviation of amplitude increases as the time increase. The dotted line 
shows the logarithm of amplitude(l-H„i„) from non-linear simulation and solid straight 

line is the same but from linear stability analysis, for ( 7 , = 0.5, w = 2.0) . 

From figure 2.7 it is clear that in the initial stages of time non-linear result of amplitude 
has good match with the linear theory but as the time increases the instability starts grow 
more rapidly than the expected value from linear theory. This is because the linear theory 

predicts the time of rupture according to [ — i.e. according to the initial potential but 

as the instability grows this term decreases rapidly resulting into smaller rupture times 
than those predicted by the linear theory. This is the reason that as the time increases the 
actual non-linear amplitude deviates from the amplitude given by linear theory. 
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Figure 2.8 Variation between the maximum value of (I), which is volume flux passing per 
unit width of the film, for ( 77 ^ = 0.5, w = 2.0) . Solid line corresponds to the linear theory 

while dotted line corresponds to the non-linear simulation. 
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0 2 4 6 8 10 


T 

Mgure 2.9. Variation between the maximum value of (v), which is the surface velocity of 
the tliin film, forr/^ = 0.5, w = 2.0. Solid line corresponds to the linear theory while dotted 
line corresponds to the non-linear simulation. 

From the above two graphs it is clear that as the time increases these values deviates from 
the linear stability theory. 

2.9.3 Effect of weissenberg number on Rupture time 

The weissenberg number has immense effect on film rupture time, in the figure below 
(Figure 2.10) it is clear, as the elasticity increases, i.e. weissenberg number increases, the 


27 




elasticity has destabilizing effect on the film stability. 



Figure 2.10. Variation of amplitude as a function of time, with weissenberg number as a 
parameter. As the value of weissenberg number increases the amplitude of instability 
increases and hence the rupture time decreases. 
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Rim height 



Figure 2. 1 1 Variation between rim height as hole diameter for ( 7 , = 0.5, w = 4.0) 

F'igure 2.10 shows the variation between rim height and hole diameter that comes out to 
be linear in nature. 
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Figure 2.12 hole diameter as a function of time, for rj^ = 0.5, vv = 4.0 

Above figure shows the variation between hole diameter as a function of time. This is 
linear in viscoelastic case. For Newtonian case, also it is linear. 

Oldroyd B model has two limiting cases: 

• Maxwell model when = 0 

• Newtonian case when t;,. = 1 , w = 4.0 
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2.9.4 Maxwell model 


1 his model is the limiting case of Oldroyd B model. In this case there is no 
solvent but only polymer is present hence the ratio of solvent viscosity to the total 
viscosity is zero. i.e. 


The constitutive relation for the Maxwell model is 


(2.53) 


4 


df. 


xy 


dt 





(2.54) 


Hence, the governing equations for the Maxwell model are same but with = 0 in the 
set of equations (2.36), (2.37) and (2.38). 

In this case also the length scale of instability is l^Jln , same as for Newtonian case. But 
time scale of rupture is decreased. 

The maximum growth coefficient for the Maxwell model is given by taking the limit of 
= 0 in the equation (growth coefficient). 


Or 


4-w 


ifw<4 


1 2(w-4) 

5- T + — lfw>4 


(w-4) 


w 


(2.55) 

(2.56) 
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Maxwell model has very interesting case when weissenberg constant approaches to 4. If 
(u'<4) then growth rate would be given by equation (above), but as weissenberg 

number approaches 4, growth coefficient increases rapidly. If (w>4) then growth 

coefficient is inversely proportional to which makes it (growth coefficient) very high. 

It is also clear from the graph below, as the value of w'= 8, the amplitude increases 
rapidly as the time increases. But for w = 2, it is not as fast as for the previous case. 



Figure 2.13. Variation of amplitude (l-/imi„) with respect to time, for different values of 
weissenberg number, for very less value of [rj^ =0.01), and comparison with the 
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LSA. Lines are from Linear Stability Analysis and points shows the results from non- 
linear simulation. 

Figure 2.13 shows the variation of amplitude with respect to time, with two different 
values of weissenberg number. One is less than 4 and other greater than 4. Ifw>4, 
growth coefficient increases rapidly for very lesser value of//,.. But it is not the case 

when w < 4 . For both these cases non-linear simulation results have shown a good match 
with the linear theory, in the initial stages of time. 

2.9.5 Purely Newtonian case 

In the Newtonian case, stress is linearly proportional to strain rate i.e. -f djV^ ^ 

where, 77 is a constant, called, viscosity of fluid. For this case there is no polymer present 
in the solution. Hence, when 77 = 77 , i.e. , 

7.=1 (2.57) 

constitutive relationship corresponds to Newto’nian fluids. 

In the constitutive equation of Oldroyd B model if we put 77 ^ = 1 then it can be simplified 
to, 

(2.58) 

So, Newtonian fluid is also one of the limiting cases of Oldroyd B model. 

In the previous work of Ghatak A., Khanna R. and Sharma A, it is already shown that 

length scale of instability for Newtonian film . And the maximum growth 

coefficient is given by, firom the equation (2.34), is independent of weissenberg number. 
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In the main governing equation (2.34) if 7 ^ is replace by one and taking weissenberg 
number to be zero its time scale is indifferent of the work done in [15] 

2.10 CONCLUSIONS AND RECOMMENDATIONS 

Elasticity plays an important role in the thin film instability'^ phenomenon. Though the 
length scale of instability remains same and hence elasticity has no influence on the 
length scale, but time scale is abruptly changes. As the elasticity or weissenberg number 
increases, time of rupture decreases and as the weissenberg number reaches its critical 
value 4 (in Maxwell Case), the growth rate increases rapidly, and hence time of rupture 
decreases. As the polymer is added in the solution, or ratio of solvent viscosity to the 
solution viscosity is decreased, rupture time is decreases, but length scale of instability is 
unaffected. Greater insight into the nature of problem can be obtained by carrying out a 
three ilimensional analysis. 



CHAPTER 3 


EFFECT OF BENDING MOMENT IN A CELL 
MEMBRANE 

3.1 INTRODUCTION 

In the process of cell adhesion, stabilizing effect of interfacial tension and bending 
moment is important. Bending moment like the interfacial tension has stabilizing effect in 
thin films. The Bending moment becomes important where there is a significant 
curvature change. In the present study, simulation work is being shown for 2-D and 3-D 
cases, for different types of fluids sandwiched between cell membrane and solid 
substrates. Both the long range Lifshitz-van-der-Waal and short range Acid-Base 
interactions are taken in the account. The non equilibrium thin films on solid substrate 
either rupture completely with the appearance of three phase contact line or undergo a 
morphological phase separation (MPS). When the apolar (LW) and polar (P) components 
of the spreading coefficient arc of opposite signs, thin non equilibrium films can undergo 
a morphological phase separation [17] Type IV films come under this category for which 
(S'-"'>() and S*“<0) 

3.2 CONCEPT OF BENDING MOMENT 

Let us consider the thin fluid film sandwiched between solid substrate and cell 
membrane. If the layer is initially flat and it is bent then a bending moment is produced 
by the expansion of upper layer and compression of lower layer, followed by an area 
change. This bending moment is related to the local curvature and expressed as [16] 


35 



M = B 


l + i-l 

V ^2 > 


(3.1) 


VVliere, B is bending rigidity constant, and is function of the membrane properties like 
compressibility of the layers. Here R^ and i?, are local radii of curvature and force due to 


bending moment is obtained asV~M . This is an extra term comes in pressure, so the 
pressure difference can be written as. 


Ap = -r 




L _L 




1 

— + - 


V^I 


R. 


2 J 


Where, y is surface tension. 


3.3 Problem Formulation 


(3.2) 


Here a fluid film is considered sandwiched between a cell membrane and solid support. 
There are two interfaces, one at the bottom of the film, between the film and the solid 
substrate and the other between the film and cell membrane at the top where film touches 
the cell membrane. 

We start with a fluid film of mean thickness ho and a small perturbation is given. The 
motion of fluid inside the film can be described by the Navier-Stokes equation, using the 
periodic boundary condition. Fluid is also considered to be incompressible in nature. 

The evolution equation, using the long wave analysis (with condition ho/L « 1) can be 
written as [13, 14], 


1 5 j ^ d / 

Updxl dx^ . 


= 0 


(3.3) 
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Where, h is the height of* free surfece, p and (j> are given by. 


2 

4 


X 



Fluid film (thickness <100nm) 


z=0 


Figure 3.1 Configuration of system 




(3.4) 


25 


IM j1 qAB 

Un 


-exp 


do-h 


(3.5) 


Where, and 5"^ are the apolar and polar component of spreading coefficient 

respectively, dp is the cut-off distance and /pis the correlation length ( = 0.6 run.) 

[sfe4.fg56] 


3.4 Non-Dimensionalisatioii 


The above equation can be non-dimensionalised using the parameters given below to 
solve for minimum number of parameters and make the equation compact. 


H = h/h,,d = djh, and / = /o//^ 


7.18 


S^'^{6d^l^) 

The non-dimensional space and time scaling, 


(3.6) 

(3.7) 
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Where 


X-ax/h^ and T=:btv/f^ 


(3.8) 


a 







V 2o^p 


o = p/ p 


(3.9) 


The equation after non-dimensionalisation can be written as 


Where , A is given by, 


(3.10) 


yK 

(6\S^^\Bd^y^ 

And, Nondimensional potential is given by, 


<I) = sgn(»S’''*') 


37JT+«exp 


d-H' 


[ I 


(3.11) 


(3.12) 


The dimensional equation is, 

(3,13) 

Equation 3.13 is the sum of four forces, per unit width of thin film, which from left to 
right are viscous forces, surface tension force, force due to bending moment and 
intermolecular force. The viscous forces merely retard the growth of instability, but 
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surface tension force tries to stabilize the thin film until the rupture occur after that it 
helps in hole growth. Bending moment term also stabilize the thin film. The 
intermolecular forces engender instability if^^ < 0 , causing the flow from thinner region 
to thicker region, thus leading to dewetting. 

3.5 Three-dimensional thin film equation 

In the three dimensional is described by introducing the y-component in the 
hydrodynamic equation. The evolution equation in three dimensional is 

(3.14) 

The three dimensional Nondimensional evolution equation is, 

li, + [//’ {A{H,, + *[h' + H„)- - H„„ - cd),.]_, = 0 

(3.15) 

3.6 LINEAR STABILITY ANALYSIS 

Linear stability analysis predicts the initial growth of thin films by linearizing the height 
around // = I . The linear equation admits solution of the form 

/i = /zo+ffexp(/fcc+(»0 (3.16) 

Where e is the ratio of initial amplitude of perturbation to mean film thickness. 

The equation can also be linearized in Nondimensional form as, 
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H = \ + Eexp(iKX + QT) (3.17) 

Where E is the initial perturbation, K is the Nondimensional wave number and Q is 
Non-dimensional growth coefficient. The linear dispersion relation relating the growth 
coefficient to the wave number can be written as, 

w = [-lit\2M)[yk^ + Bk^ + (3. 1 8) 

It can also be expressed in the Nondimensional form as 

^ = + +(<!>« X,.,] (3.19) 

i.e. 

Q = -K'^-sgniS^''') 1-Pexp > K' (3.20) 


Growth coefficient will increase as the wave number is increased. It reaches a maximum 
and then starts decreasing, so the dominant wave number is defined as km, for which 

growth coefficient co is maximum, or — = 0 . 

dk . 



The critical wavelength is defined the wave number at which growth coefficient vanishes 
i.e. £y = 0 . 
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i.e. for the wave length greater than the critical wavelength, all perturbation will die out 
and h = will the steady state solution. 

In the Nondimensional form it can be written as, 



(3.23) 



2 


(3.24) 


There are two asymptotic cases. 
1)5 — > 0 


-mr 


-(A), 




ly 


1/2 


(3.25) 


2) y -^0 


^mB ~ 


■(A) 


A-Ao 


35 


1/4 


(3.26) 


The dominant wave length is related to the dominant wave number by the relation. 


A=2n/K^ 


(3.27) 


The time of rupture can also be predicted by the linear theory, by putting H 0 and 
cos(KX) = -1 in equation (3.17). This gives the time of rupture, i.e. the time at which first 
time film surface touches the solid substrate, in terms of initial perturbation amplitude 
and maximum growth coefficient. 
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(3.28) 


T = 



_ 1 _ 

E 


3.7 Numerical Method (2D) 

The non-linear equations, derived earlier are sixth order nonlinear partial differential 
equation in space. It is discretized using central differential technique using half node 
interpolation. The resulting set of ODEs are very stiff in nature because the wide 
disparity between length scale and time scale. Hence this initial value problem is solved 
using Gear's technique which is appropriate for stiff set of differential equations. Nag 
library subroutine D02EJF which employs Gear’s algorithm is used to integrate the ODEs 
Equation (3.10) is solved in an X domain using the periodic boundary condition. This 
length is divided! in to N parts using the finite difference hence this PDE is being 
converted in to N number of ODEs, each node represents an ODE. 

A.\' = L/(7V-1) = A 
Wc can reach at /'*' iukIc by 

A'=(/-1)A (3.29) 

If the film thickness at the node is H, then derivative at this point is given by 


dH, 

dT A 


(3.30) 


Where, 

f 1, = [ H + ®),v ], = (3.31) 

F2, +.®], 
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Where, this spatial derivatives are further discretized using the central difference 
technique, given below. 


[«.»■],= 




(3.32) 


And 


[H,v,,xvl 




(3.33) 


3.8 Numerical Method (3D) 

The three dimensional equation is solved in domain using the periodic boundary 

condition. The equation is discretized using the central difference technique using the 

half node interpolation. Because the equation is very stiff in nature hence D02NCF 

2 

routine from the Nag library, which employs Gear’s algorithm, is used to solve the N 
number of differential equation. Fluid surface is parallel to the X-Y plane and Z direction 
is normal to the film. Finite difference technique, taking both the difference in X and Y 
direction same, has been used. 

AY = Ar = I/(Y-l) 

To reach at any location at (Xi, Yj) we can simply write 

(Jf„J',)=((/-l)A,p-l)A) (3.34) 

If the film thickness at the (Xj, Yj) location is Hij then the derivative at this point can be 
written as, 

^ y — 2 2 3 N (3 35) 

dT ~ A A 

Where function FIX is given by, 
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/•'l.V 


// •( A ( ^ +//),.)+ + H + ^>} = H] 

•' -ii.i 


F2 -F2 

3 •' •“ 


A 


(3.36) 


Similarly the function FIY is given by. 


f\Y 


XVXV ^YYYY ^}yl. . 

J/,/ 




FI -F2 


iJ-\/2 


A 


(3.37) 


Where F2sj is given by, 


I'l., ^ (3.38) 

Discretization in space is done, using cenntral difference technique in space, the whole 
domain is devided in to N x N number of points. Now at each point the partial differential 
equation converts to an ordinary difference eauation. hence this PDE in the whole domain 
converts to the N" number of ODEs. Nag library subroutine D02NCF is used to solve 
these N^ number of ordinay differential equations. This subroutine employs Gear’s 
technique which is best souted for the stiff differential equations. The 60 X 60 points 
were found to be sufficient for the grid convergence test. 


3.9 RESULTS AND DISCUSSION 

The total spreading coefficient can be calculated using the equations below for LW and 
Acid-Base components, for the system described above. 
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^AB 


= 9 



fc -^|K 




(3.40) 


\Vheie, subscript c, t, s are corresponding to the cell membrane, fluid film and the 
substrate. S and S are spreading coefficients for Lifshitz— vander Waals and Acid- 
Base components, respectively. The typical value of bending moment is taken to be 10’' '^ 
ergs [7] for all the cases. 

The relative importance of bending moment is clear if the dominant length scale of the 
case bending moment with no interfacial tension is compared with only interfacial 
tension case. 

The dominant length scale for the case with bending moment and no interfacial tension is 
given by, 




3B 




h^ho 


1/4 


(3.41) 


Similarly, the dominant length scale for the case with interfacial tension with no bending 
moment is given by, 


^mr ~ 


2y 


-(aL, 


1/2 


The ratio can be written as. 


/'t d/ j. \ 




3SK), 




4f 


(3.42) 


(3.43) 


It is clear from equation (3.26) that for low values of surface tension, the effect of 
bending moment has to be considered, otherwise length scale of instability will be very 
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low and hence long wave analysis will not be applicable in this case. Lastly, for higher 
value of ho and y, the effect of bending moment can be neglected eq. (3.18). The bending 
moment stabilizes the system as it is evident from figure (3.2) which shows the inverse 
relation of growth coefficient with the bending moment. 



fig. 3.2 Variation of growth coefficient with the bending moment. 


To analyze the system, it has been dividend in to two main parts. First is the case when 
the Acid-Base interaction vanishes, i.e. = 0) and hence (P = 0). The second part is 
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when ( P ^ 0 ), this part is further subdivided into two parts. The first sub-case is when 
the substrate is apolar and in other sub-case substrate is polar. 

CASE 1: P = 0 

In this case, the fluid film and any one of, either the substrate or the cell membrane is 
apolar in nature. For this case the Acid-Base component of the Spreading coefficient 
vanishes, i.e. S '''* = 0, it is also self exemplified from the equation (3.40). Hydrocarbon 
film in contact with the lymphocyte cell is one of the examples. component should 
be taken to be negative to describe the attractive behavior of Lifshitz-vander Wall 
component. 

For this caste from the equations (3.12) and (3.23), it is clear that dominamt wave number 
has inverse dependence on A, a function of y and ho- The behavior for this case is of Type 
1 system, where is negative for all film thicknesses. It is observed from the non-linear 

simulations that for any value of B, the system will be destabilized from the effect of 
attracting (destabilizing) Lifshitz-vander Wall component [10]. 

Case: 2 (P ^ 0) 

In this case, both the polar and apolar interactions are taken into account. This system is 
further divided into two parts, in the first one apolar substrate is taken and in the second 
one polar substrate is taken. In both the cases, cell membrane is taken to be polar in 
nature. In the first case, when the substrate is apolar, the system is chosen in such a way 
that component is positive i.e. fi:om the equation (3.39) it is clear that for to be 
positive (stabilizing), the Lifshitz-vander Wall component of the free energy of the film 
should lie in between the LW component of free energy of substrate and the cell 
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membrane, i.e. y^. > y‘f > y^'^ or y‘f <y‘f' Kyf , and component of free energy 

to be negative (destabilizing). In the second case substrate is chosen in such a way that 
component is negative and component is positive. The above two cases comes in 
the category of Type II system and Type IV systems. 

Type II Systems 

Type II systems are a combination of long wave repulsion and short wave attraction. 
which is long range force, is negative while S^, which is short range Acid-Base force, is 
positive. 

Type IV Systems 

Type IV systems are a combination of negative, long range attraction of Lifshitz-vander 
Wall forces and positive, short range repulsion of Acid-Base forces. The graphs for these 
two systems are shown below. 
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2 4 


Fig 3.3. Variation of wit: 
mJ/m" and S'''** is -33.0 mJ/i 



Fig 3.4. Variation wi 
mJ/nr and = 25.0 mJ/r 




Apolar Substrate 

Figure 3.5 shows the profile of thin film in 2-D case, in 2A length scale where A is the 
dominant wave length corresponding to that film thickness, with mean thickness as the 
parameter. As the initial mean film thickness is increased, there is the increase in the peak 
value. The Nondimensional time of rupture is increased as the mean thickness is 
increased, which indicate that film rupture is faster in thinner films. 

Figure 3.6 shows the effect of interfacial tension at constant mean film thickness. It 
shows that as the value of interfacial tension is increased there is an increase in the 
rupture time; this is because surface tension has stabilizing influence. 

Figure 3.7 shows the whole morphology of thin sinusoidal film, with mean thickness ho = 
5 nm and initial perturbation of e, on a length scale of 2A, where A is the dominant 
length scale. Because of the negative value of which signifies, higher pressure in 
thinner region and lower pressure in thicker region. Hence the flow starts from thinner 
region to thicker region leading to the film stability. The amplitude of wave start growing 
and at one point the film surface touches the surface and making hole which starts 
growing with time. Simultaneously a rim starts growing near the hole, which also 
increases as the time increases. Finally two neighboring rims start interacting and make a 
bigger drop. 

In Fig 3.8, length scale of 3A is taken for the simulation, where A is the dominant length 
scale given by linear stability analysis. A thin film of mean thickness of 5.0 nm and initial 
amplitude of random perturbation of (8=0.01) was taken for the simulation. Three waves 
start forming on this length scale which validates the simulation. 
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Figure 3.5. Final profile of the thin film, 
sinusoidal perturbation, for apolar subsl 
Surface tension y = 0.001 mJ/m^. 
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Figure 3.6 Final profile of the thin film, with surface tension as the parameter, starting 
with sinusoidal perturbation, for apolar substrate, = 2.0 mJ/m^ and = -33.0 
mJ/m“. Mean height ho = 5 nm. 
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Figure 3.7 Full morphology of thin film starting with mean film thickness ho = 5 nm. On 
2 A length scale, for apolar substrate = 2.0 mJ/m^, S'*"® = -33.0 mJ/m'. Surface 
tension y = 0.00 1 mJ/m^. 



Fig 3.8. Morphology of thin film starting with mean film thickness ho = 5 nm and with 
random perturbation of s = 0.01. On 3A length scale, for apolar substrate = 2.0 
mjW, = -33.0 mjW. Surface tension y = 0.001 mJ/m^. 
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Polar Substrate 


In this case subsltaic is taken to be apolar in nature. The values of spreading coefficients 
are taken such that this s>stem behaves like a Type II system i.e. S‘-'^ is negative and 
is positive. Again the sign of both the spreading coefficients are negative, we expect 
Morphological I’hase Separation (MPS) i.e. formation of a stable interconnected 
microstructure consisting ot an array of small drops in equilibrium with the thin films. 
Figure 3.9 shows the final profiles of thin films on polar substrate on 2A length sc ale 
starting with a sinusoidal initial perturbation, taking mean height as the parameter. As the 
mean height is increased, rupture time as well as the peak height increases. 



Figure 3.9. Profiles of thin films on an apolar substrate, varying mean height as the 
parameter. = -4.0 mJ/m^ and = 25.0mJ/m^. 


54 



Figure 3.10 shows the 3-D simulation results for 5 nm thick film on a homogenous 

substrate. The length scale of instability is recognized by the linear stability analysis. It 

was found that two different morphological patterns and their sequence of evolution by 

which dewetting can occur, depends on the form of the potential in the neighborhood of 

the initial film thickness. Figure 3.11 summarizes the major events in the time evolution 

of pattern in a relatively thick 80 nm film, which is far more than the height at which 

(j),^ shows minimum. The initial random disturbance is first reorganized into a small 

} 

amplitude bi-continuous pattern on a length scale close to A . Long “hills” of the structure 
undergo some fragmentation, while the “valleys” thin locally to produce largely circular 
holes surrounded by circular uneven (in height) rims. The general pathway of evolution, 
whenever the initial mean thickness is sufficiently high so that the repulsive interactions 
at the minimum thickness are encountered only after a considerable growth of instability, 
can be summarized as (bi-continuous pattern circular holes with rims hole 
expansion -> hole coalescence) [8]. At the later stages of time it shows the breaking of 
rims and formation of some ridges like pattern, this instability is called Rayleigh 
I Instability. The dewetting occurs by expansion of holes. 

In contrast to the above scenario, figure 3.12 depicts a different pathway of evolution for 
the relatively thin films which encounter repulsion in the early stages of growth. A bi- 
continuous structure composed of long hills and valleys persists until the fragmentation 
of “hills” producing an array of microdroplets. These droplets become increasingly 
circular due to surface tension and increase in height due to flow from “valleys” which 
thins and flatten out. The droplets are circular due to the surface tension force. In this 
case, the dewetting occurs by the retraction of droplets, rather than by expansion of holes. 
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Figure 3.10 Different stages of evolution in 5.0 nm thick film, starting with random 
perturbation, lo = 0.6 nm, SLW = 2 mJ/m^ and = -33.0 mJ/m^ the area of each box is 
1 6A^ (A = 2.98, is the dominant length scale) the images correspond to from left to right 
at time T = 0, 0.0016, 0.0060, 0.0130, 0.0180, 0.0215, 0.0225, 0.0228, 0.0236, 0.0240, 
0.0248, 0.0277 respectively. Increasing shades indicates increase in height. 
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Figure 3.11. Different stages of evolution in 80 nm thick film, starting with random 
perturbation, lo = 1.0 nm, SLW = -4 mJ/m^ and = 25.0 mjW the area of each box is 
16A^ .The images correspond to from left to right at time T = 0, 1.08, 14.32, 53.39, 
59.5 1 , 63.02, 66.06, 7 1.86, 76.99, 95.20 respectively. Increasing shades indicates increase 
in height. 
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Figure 3.12 Different stages of evolution in 3.5 nm thick film, starting with random 
perturbation, lo = 2.5 nm, SLW = 4 mJ/m^ and S'^*^ = -0.1 mJ/m^ the area of each box is 
16A^ the images correspond to from left to right at time T = 0, 128.46, 472.07, 587.67, 
629.12, 718.35, 940.31, 1811.97 respectively. Increasing shades indicates increase in 
height. 

3.10 CONCLUSIONS AND RECOMMENDATIONS 

Bending moment plays an important and dominant role in the thin film instability. It 
plays an important role when surface tension has very low value. Length scale and time 
scale is changed, while the time scale has not been affected much by the inclusion of the 
bending moment because this instability is a micro drop phenomenon, but length scale is 
changed largely. The 3-D figures show that different initial height shows the different 
patterns and different phenomenon of dewetting, depending upon after what depth it feels 
the repulsive forces. 
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